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Abstract 

We quantum mechanically analyze the fractional quantum Hall effect in graphene. This will 
be done by building the corresponding states in terms of a potential governing the interactions 
and discussing other issues. More precisely, we consider a system of particles in the presence of an 
external magnetic field and take into account of a specific interaction that captures the basic features 
of the Laughlin series v = 2fp[- We show that how its Laughlin potential can be generalized to 
deal with the composite fermions in graphene. To give a concrete example, we consider the SU(N) 
wavefunctions and give a realization of the composite fcrmion filling factor. All these results will 
be obtained by generalizing the mapping between the Pauli-Schrodinger and Dirac Hamiltonian's 
to the interacting particle case. Meantime by making use of a gauge transformation, we establish a 
relation between the free and interacting Dirac operators. This shows that the involved interaction 
can actually be generated from a singular gauge transformation. 
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1 Introduction 



Nowadays the observation of the famous quantum Hall effect (QHE) [lj does not remain at the 
stage of the semiconductors, but it can be seen in other materials. This is for example the case 
of graphene, which is a projected graphite of the group symmetry C3 into two-dimensional spaces. 
When such system is submitted to a perpendicular magnetic field, it appears an anomalous integer 
Hall conductivity [2j [3] that is a manifestation of the relativistic particle motions. It is originated 
from different sources, which are exactly the four- fold spin/valley and the Berry phase due to the 
pseudospin (or valley) precision when a massless (chiral) Dirac particle exercises cyclotron motion. 
This result was theoretically suggested by two groups [U [5] independently. 

The emergence of QHE in graphene during 2005 offered a laboratory for different investigations 
and applications of some mathematical formalism. Many of them are brought from the QHE studies 
in semiconductors and used the Dirac operator as the cornstone instead of the Landau Hamiltonian 
as the starting point. Indeed, different questions have been partially or completely solved but still 
other attentions to be paid. This concerns an eventual fractional quantum Hall effect (FQHE) [6j in 
graphene and related issues. Such interest can be linked to the fact that FQHE is a fascinating subject 
because starting from its appearance is still capturing a great attention. This is due to its relations 
to different areas of physics and mathematics. 

There are many interesting theories have been appeared dealing with some problems related to 
the anomalous FQHE, for instance one may refer to the papers [10]. In particular, based 
on different arguments some wavefunctions have been suggested as good candidates to describe the 
subject [TTJ [T21 [T3] . For instance, some of them used the composite fermion approach [T3] and other 
made a straightforward generalization of the Halperin theory [15] for non-polarized spin. However, 
some of them can be linked to each other and recovered from some general proposals as has been seen 
in [13] as we will see in the present work. 

Another theory was proposed by the first author [T3] who studied the effect exhibited by the 
relativistic particles living on two-sphere § 2 and submitted to a magnetic monopole. In fact, he 
started by establishing a direct connection between the Dirac and Landau operators through the 
Pauli-Schrddinger Hamiltonian Hf p . This was helpful in the sense that the Dirac eigenvalues and 
eigenfunctions were easily derived. In analyzing Hf p spectrum, he showed that there is a composite 
fermion nature supported by the presence of two effective magnetic fields. For the lowest Landau level 
(LLL), he argued that the basic physics of graphene is similar to that of two-dimensional electron 
gas, which is in agreement with the planar limit. For the higher Landau levels, he proposed a SU(N) 
wavefunction for different filling factors that captures all symmetries. Focusing on the graphene case, 
i.e. N = 4, he gave different configurations those allowed to recover some known results. 

Motivated by different analysis concerning the anomalous FQHE and in particular [T3], we de- 
velop a real approach in order to describe the subject. Taking advantage of our knowledge about 
the phenomena in semiconductors, we introduce some effective interactions as relevant ingredient in 
forming the Hamiltonian system. For this, we consider M-particles in the presence of a perpendicular 
magnetic field B where the interaction is taken to be two and three-body types. These have been suc- 
cessfully applied to the conventional FQHE, for instance one may see |16tll7]. Based on the connection 
between the Pauli-Schrddinger and Dirac Hamiltonian's, we construct an appropriate Dirac operator 
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that captures the basic feature of the interacting system and where the Laughlin wavefunction [18J 
is its groundstate. This has been used to derive a basis as set of the excited states as well as the 
corresponding energies. 

Because of the basic physics is the same in both systems: semiconductors and graphene for 
the Laughlin series, we start by establishing the corresponding effective potential. It will be gen- 
eralized to describe the composite fermions in graphene [19] of filling factor u c f = 4 (re + 5), with 
n = 0, ±1, ±2, This of course captures the physics behind the fractional filling factors beyond 
Laughlin series 2j+T> with I is integer, for FQHE in graphene. Moreover, we return to the SU(N) 
wavefunctions [13] to firstly give for the step a concrete example showing a particular type of poten- 
tial. Secondly, we show how the corresponding filling factor can be linked to that generated from the 
composite fermions pictures |19j . Different discusions will be reported elsewhere about the matter. 

The present paper is organized as follows. In section 2, to show the difference between the in- 
teracting and free particles cases, we start by studying a system without interaction. We use an 
approach based on the mapping between different involved Hamiltonian's to get the corresponding 
spectrum and eigenfunctions. In section 3, for the later convenience, we give some discussions about 
the Laughlin potential that captures the interaction effect and is behind the filling factor v = ^jxj • We 
introduce some specific interactions that are the two and three-body type and analyze the behaviour 
the new system in section 4. In particular, we show that the Laughlin wavefunction is a groundstate 
of the present system and resort the spectrum's for the excited states for the three Hamiltonian's. 
The considered interactions can be generated from a singular transformation and allow us to establish 
a relation between both Dirac systems: free and interacting. In section 5, we consider the composite 
Dirac fermions in terms of an appropriate Jain potential, which can be obtained by making use of 
straightforward generalization of the Laughlin one. In this framework we analyze the SU(N) wave- 
functions and show its basic features in section 6. Finally, we conclude and give some perspectives of 
the present work. 



2 Dirac Hamiltonian 



We start by developing our approach that will be subsequently generalized to the case of the interacting 
system. This concerns a mapping between different spectrum's: Landau, Pauli-Schrodinger and Dirac. 
Indeed, let us consider one-relativistic particle living on the plane (x, y) in presence of a perpendicular 
magnetic field B. This can be described by 

1 r / _ e _\ -1 2 

H P s = ^~ a-(P--A) (1) 

2m 1 \ c /J 

where a are the Pauli matrices and verify the usual relations, namely 

{crj,<7j} = 25ij, [<Ti,<Ti] = 2eijk(Jk- (2) 

We will see how the Pauli-Schrodinger Hamiltonian Hps can be used to get the Dirac operator for the 
present system. Specifically, it will be derived from the square root of -Hps- This mapping is useful in 
sense that the Dirac eigenvalues and eigenfunctions will be easily obtained. 

One way to establish the mentioned mapping is to express Hps in terms of the Landau Hamiltonian, 
which describe free particle. This connection has an important advantage because it leads to exactly 
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derive different quantities related to the subject. Thus, by choosing the symmetric gauge 

A = ^(-y,x) (3) 
we show that Hp$ can be mapped in the form 





Hps = : TT - I I (4) 



where Hi is simply the Landau Hamiltonian and reads as 

H L = — (P- -AY (5) 

2m V c / W 

which is extremely used in different contexts and in particular in analyzing QHE in the semiconductor 
systems. To easiest way to resort different spectrum is to write Hi in terms of the annihilation and 
creation operators as 

ff L = 1(0*0 + 00*) (6) 
where, in the complex notation z = x + iy, a and a) are given by 

a = 2 ! + f* a, =- 2 i+f- < 7 > 

They verify the commutation relation 

[a,a Jf ] = 2B. (8) 

Since we are wondering to generalize different Hamiltonian's entering in the game, let us convert Hi 
into its analytical form. It is not hard to obtain 

Hereafter, we set the fundamental constants (e, c, H, m) to one. Clearly, the spectrum and eigenfunc- 
tions of Hps can be derived from that of Hi. 

The above tools can be applied to analysis the anomalous QHE in graphene. Indeed in such 
systems, the two-Fermi points, each with a two-fold band degeneracy, can be described by a low- 
energy continuum approximation with a four-component envelope wavefunction whose components 
are labelled by a Fermi-point pseudospin = ±1 and a sublattice forming an honeycomb. Specifically, 
the Hamiltonian for one-pseudospin component can be obtained from ([4]) under some consideration. 
This is 



i 



a) 



H " = T2 Vf \-a )■ (10) 
where Vf ~ is the Fermi velocity, which will be set to one, and the many-body effects are neglected. 
Its spectrum can be determined in a simple way if we introduce its square. This is 

g = 1 / a*a y 

D 2 ^ aa) J V ' 

It is related to the Pauli-Schrbdinger Hamiltonian up to some multiplicative constants. It is clear 
that, Hq is written in terms of the diagonal form of the Landau Hamiltonian ([9]). Therefore, it 
spectrum can easily be obtained. 
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We start by solving the eigenvalue equation 

H&V = Efa (12) 

Since has to do with the Landau Hamiltonian ([9]), the wavefunctions \E' should be written in an 
appropriate form. This is 

*m,n = ( ^ ) (13) 



m—l,n 



where the eigenfunctions ip m n are given by 



^ m , n (z, z) = ) ' - = z n L^( T )e~T»*, m, n = 0, 1, 2 ■ ■ ■ . (14) 
y/2 n+1 ir(m + n)l ^ 

The corresponding Landau levels are given by 

El(m) = Bm. (15) 
We can show that the normalized eigenfunctions of Hq take the form 

*Wo,» = 4( * / H ' n ) (16) 

V2 \ -lip\m\-X, n J 
with the convention ifr—i n = 0. Note that, the zero mode wavefunction is 

y(0,n) = ( V-0,n \ _ (17) 







Their energy levels read as 



E D (m) = ±y/B\m\. (18) 



This makes a difference with respect to the Landau spectrum for the present system. In fact, it has no 
zero energy and is discrete as well, unlike (I18p . These what make the integer QHE is an unconventional 
effect in graphene. 

The above analysis can easily be generalized to many-body system without interaction. This leads 
to a spectrum as sum over all single particle ones and eigenfunctions as tensor products. However, we 
will give a generalization of the system by introducing a kind of interaction that is behind FQHE in 
semiconductors and see what is its influence on graphene. 



3 Laughlin wavefunction 

To develop our main idea we need first to start from what we know so far about the Laughlin wave- 
function. This latter is involving a Jastrow factor that is originated from a specific interaction between 
particles. In fact, it corresponds to an artificial model that captures the basic physics of the filling 
factor v = 2fTi- This is an typical example to get more general results, which concern the composite 
fermions as well as the SU(4) wavefunctions for graphene. For our task it is necessary to start from 
the Laughlin theory for FQHE in graphene. 
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To talk about the Laughlin series for FQHE, it is convenient to introduce some discussions about 
LLL because it is a rich level and has many interesting properties. Indeed, in such level the Landau 
system behaves like a non-commutative one, which is governed by the commutation relation 

M]lll = §. (19) 

This effect can be interpreted as follows. In the present level, the potential energy is strong enough 
than kinetic energy and therefore the particles are glue in the fundamental level. Using the same 
analysis as for the case of H^, is not hard to get a basis as set of the eigenstates. They are 



1^7 LLL = 

V m' 



= (a t |LL L r|0)LLL. (20) 



where the annihilation and creation operators reduce now to 



2 ' 1 2 

The projection into the complex plane leads to the eigenfunctions 



B B 
o|lll = — z, a) | LLL = —z. (21) 



$ / z g\ = .j " Z ™ e -#M 2 . (22) 

v ' ' V 2 m+1 vrm! v ' 

As we will see very soon, this can be generalized to get more interesting results. In particular, these 
will offer a way to get in contact with the Laughlin theory [18] for FQHE. 

Returning now to discuss FQHE in graphene. In doing so, Let us consider M-particles in LLL, 
which of course means that all = with % = 1, ■ ■ ■ , M and each mj corresponds to the spectrum f) 161 
[T8|) . The total wavefunction of zero-energy Landau level (|17p can be written in terms of the Slater 
determinant. This is 

rftz, -z) = j^z™ ■ ■ ■ z^ 1 exp U| f) \ Zi \^j (23) 

where e n '"* M is the fully antisymmetric tensor and rrii are integers. It is relevant to write this wave- 
function as Vandermonde determinant. We have 

M / B M \ 

ip(z, z) = const (zi - zj) exp ( -— ^ \zi\ 2 J . (24) 

i<j \ i / 

The corresponding energy is 

E(rm) = B m i + ■ (25) 
This can be interpreted by remembering the Laughlin wavefunction 

AI / M \ 

4augh(^, z)=H (zi - Zj ) 2l+1 exp -- N 2 • (26) 



i<j 



It is well-known that it has many interesting features and good ansatz to describe the fractional QHE 
at the filling factor v = 2iqrr> with I is integer value. It is clear that (|24p is nothing but the first 
Laughlin state that corresponds to v = 1. Actually, (I2"4l is describing the first quantized Hall plateau 
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of the integer QHE. It is also convenient for us to write (|26p under another form. This is Note that 
can also be written as 

IOLau g h = { e i ™C---^ M r +1 |°)" ( 27 ) 
Consequently, the wavefunction for the particles in the graphene LLL are identical to those of the 
2DEG LLL. We conclude that the basic physics in both systems is the same. More discussion about 
this issue can be found in [9]. 

As far as we know FQHE is a consequence of the interacting particles in the presence of an external 
magnetic field. Then, obviously the Laughlin wavefunction is a manifestation of a specific interaction 
type, which coded in the Jastrow factor. At this point, one may ask how does look like such interaction 
and to reply this question, let us write (f26|) as 

iLrtfoz) =F AM *{z,E)exp i-j J>| 2 ) (28) 

where the function F Laugh (z, z) ant its anti-holomorphic partner are given by 

M M 
F Laugh (z, z) = H fa - z 3 ) 2l+l , F Laugh (z, z)=H (zi - -z 3 f l+1 . (29) 

These are nothing but the results of a specific interaction, which generated from particles. To clarify 
this point, let us define two functions as the derivatives with respect to the variables z and z of F Laugh 
and F Laugh , respectively. They are 

(V; Laugh ) < = -2i^-(lnF Laugh - lnF Laugh ), (Vt augh )i = 2i^-(lnF Laugh - lnF Laugh ). (30) 

Explicitely, the corresponding potential can be written as 

(yUugh ) . = _ 2i(2Z+1) ^_L_ j (^Laugh )j = _ 2i(2Z + 1) J-_^_ (31) 

Zi Z i Zi Z n 

which is resulting from of interaction of the i th particle with the remaining j ones. These can also 
be identified to the components of the Chern-Simon potential, which has been employed at many 
occasions to study the conventional QHE. For instance, more discussion about such identification can 
be found in [T7j. This of course open another way to describe the anomalous QHE by making use the 
field theory approach. Strictly speaking the present potential is, in reality, a part of an Hamiltonian 
whose groundstate is exactly the Laughlin wavefunction. Moreover, we show that how (|30p can be 
generalized and used to discuss other filling factors beyond Laughlin states for FQHE in graphene. 
This mainly we will be done by adopting the composite fermion picture [13] and using the SU(4) 
wavefunctions [13] as an illustration. 

4 Many-body interactions 

As claimed before we are wondering to investigate the basic feature of FQHE in graphene. In doing 
so, we have to take account of a kind of interaction. For this, we start with a Hamiltonian that has 
been considered at different occasions and is relevant to describe FQHE at the filling factor v = 
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Further, we show that this will be applied to graphene and allow us to derive different results. Mainly, 
we apply the approach seen before for free relativistic particles to achieve our goal. In fact, we start 
from a Hamiltonian generalizing the Landau one to map a bridge between the Pauli-Schrodinger 
and Dirac Hamiltonian's for interacting particles. Of course, this will be done by constructing the 
appropriate annihilation and creation operators those will be used to build the excited states and 
derive the corresponding energies. 

4.1 Excited states 

In the beginning, we start by resorting the spectrum for a system that include a specific interaction. 
This can be done by constructing an explicit model that has the Laughlin wavefunction as an exact 
groundstate. Let us particles living on the plane in the presence of an external magnetic field and 
involve some kind of interactions. When these are of two-body and three-body types, the present 
system can be described by the Hamiltonian 



i M M i / a \ M i / a 

\ E <* " M? + V E TT^TT ~ ~ V E T^TT (4 " «* 



(zi - Zj) \ dzi J fri. (zi - zj) \ dzi 



where the coordinate Zi = x% + iyi denote the position of the i th particle and r\ is the inverse of the 
Laughlin series, i.e. r/ = 21 + 1 is an odd integer. The second and third terms in (|32p are two-body 
interactions whereas the fourth one is a three-body type. In the symmetric gauge where the complex 
components of the vector potential A are 

A z = - % -Bz, A- z = % -Bz (33) 

the Hamiltonian (1321) takes the form 



M < d 2 „ / a a\ b 2 l i ( d b 



Hr = 4^{ 4 ^- s (^-^)-T^ + ^(^)l 2 ^"^ 

A 1 / d B \ 2 ^ 1 

It is obvious that when 77 is switched off we get the Landau Hamiltonian for a free system of M- 
particles. Therefore, the presence of this parameter will play an important role in discussing different 
issues related to the anomalous FQHE. 

To analyze the behaviour of interacting relativistic particles, it is relevant to derive the correspond- 
ing Dirac Hamiltonian. As far as we know, the easiest way to do is to adopt the technical used in 
section 2, which based on the mappings between three Hamiltonian's. More specifically, this can be 
done by making use of an algebraic method that requires to introduce the appropriate annihilation 
and creation operators. Indeed, in terms of the interaction they can be realized as 

i d B ^ 1 it d B ^1 

CJ Zj , Zi Zq Cry ■ ^ , Zj Z 1 
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It is not hard to show the commutation relation 

M M 



2B6ij - 4rj S{Zi - Zj ) - 4r? J] 5{ Zi - z 3 ) (36) 



where for instance the second term is resulted by acting d Zi on the third term of and so on. These 
allow us to get a diagonal form of H'" 1 , such as 



rrint 

2 



I f M \ M 

- 4 A + MB\-2r 1 Y J 5( Zl - Zj ). (37) 

\ i / i+j 



Since 6(zi — zj) is irrelevant when acting on an antisymmetric wavefunction, it can be ignored. More- 
over, Ai and M satisfy the canonical harmonic oscillator commutation rules even for non-zero rj. 
Hence, the Hamiltonian is exactly soluble for all the excited states. A basis of wavefunctions describ- 
ing excitations of the v = ^ states in the Laughlin theory is therefore given by 

\{m i }, V ) L =e^A{ 1 ...Al M .^t ) |0) (38) 



where are integers such that < mi < • • • < tum ■ The corresponding eigenvalues are 

4 nt (W,7?) = B^m i + ^M[(7 ? -l)(M-l) + l]j . (39) 

For later convenience, let us project these states on the complex plane. Indeed, by requiring that 
m = i — 1, we end up with the wavefunction 

M M / R M \ 

4 (z, z) = U(z k - Zl ) n ( Zi - z 3 ) 2l+1 exp -- l^l 2 • (40) 

k<l i<j \ i / 

An interesting result can immediately be noticed. Indeed, now it is obvious to see that the Laughlin 
wavefunction is an exact groundstate of ([M]) and corresponds to the energy H^-. More precisely, it is 
easy to check that Ai annihilates Laugh = 10 Laugh) namely 

A|»7) Laugh =0, Vi. (41) 

This is nothing but the lowest Landau level condition. It is clear that for 77 = 1 we recover the fully 
occupied state seen before. This show that the present model captures some basic features of the 
phenomenon. 

In the next, we will show how the above tools can be used to describe FQHE in graphene. In 
the beginning we start by establishing the appropriate background to do our task. In particular, we 
determine the Dirac Hamiltonian and related matters for the present system. 

4.2 Spinors 

To apply the technical used before, it is convenient to write down the Pauli-Schrodinger Hamiltonian 
-ffpg for interacting particles in the presence of B. The easiest way to do is to give its expression 
in terms of the Hamiltonian H'" 1 . The advantage of this mapping is make use an algebraic method 
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to diagonalize and therefore derive the associate spectrum. This will be helpful in establishing the 
corresponding Dirac formalisms. 

One way to write down Hp£ is to generalize the Hamiltonian for free system to the interaction 
case. In doing so, we get the form 

This is actually involving three matrices and the second one is particular. On the other hand, it can 
be identified to that used by Koshino and Ando [21] in discussing the diamagnetism in disordered 
graphene in arbitrary magnetic fields. However, in our case it be ignored since does not affect the 
antisymmetric wavefunction. 

To start talking about different issues, we need first to resort the spectrum as well as the eigenstates 
of Hp$. Clearly, these matters can be easily derived to get the excited states as spinors 

IKU> PS =( limU)L ) (43) 
V \{mi -1},v)l J 

where the eigenstates \{rrii},rj)i are those for the Hamiltonian (|37p . with the convention |{— l},r/) is 
a null state. The corresponding eigenvalues are 

E&(m, V ) = B^m i + ^(r,-l)M(M-l) S J . (44) 

We have some remarks in order. For mj = and 77 = 1, there is a zero mode energy that corresponds 
to a fundamental state. However, for the case where m, = and 77 7^ 1, we end up with 

E^(m = 0,7]) = -(17 - 1)BM(M - 1). (45) 

This captures the basic features of the excited states for the present system, which is exactly the 
LLL condition. It offer another discussion about spin as a degree of freedom and its influence on the 
discussion reported before. Specifically, its relation with the conventional QHE, but this is out of 
scope of the present work. 

At this level more discussion can be reported on the spectrum and in particular the corresponding 
Hilbert space, which is a set of all states \{mi},r]}ps- It is invariant under the semi-direct group 
G = U(l) x C and therefore for an unitary projective representation T(g) of G we have the invariance 
property 

T(g)Hl"i = H^T(g), \/g G G (46) 
where an element of G is a 2 x 2 matrix of the form 

«=(;;)■ (47) 

This representation can be further explored and then make a group theory approach to investigate 
the basic features of the present system. It is acting on the spinors trough the relation 

T(g)\{m}, V ) PS =j(g,z)g- 1 \{m},T 1 ) PS (48) 
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where the quantity g 1 |{w-},??)ps is the pull-back of |{m},r/)ps through the mapping z \— >• g ■ z such 
as 

g ■ z = az + b. (49) 

The function j(g,z) is and automorphic factor. One also can add more analysis about the incom- 
pressibility in terms of the area preserving diffeomorphism algebra that can be realized in terms of 
the ingredients of the present system. 

4.3 Interacting Dirac particles 

The above mathematical background can be applied to analyze FQHE in graphene. For this, it is word 
while to establish a Dirac formalism for the interacting particles. One way to do this is to make an 
analysis based on the former approach for non- inter acting particles. This requires an analogue Dirac 
operator in terms of the annihilation and creation operators. 

Right now we have all materials needed to do our job. Indeed, the appropriate interacting operator 
can be written in similar way to (|50p . Thus for the i th particle, we have 

- in ( -\ ? ) (50) 

where A{ and A- are given in (|35p . To make a comparison with the Pauli-Schrodinger Hamiltonian, 
let us define the square of HS 1 . This is given by 

Summing over all M-particles and if we forget about all constants entering in different Hamiltonian's, 
it is not hard to end up with the relation 

(#b nt ) 2 = H$£. (52) 

This equivalence shows that how the applied approach is relevant as well to analyze the interacting 
particles. Remember that 8{zi — zj) and 5(z~i — z~j) can be dropped because they do not affect the 
antisymmetric wavefunction. Therefore, we can carry out our previous study concerning free Dirac 
particles and its connection to FQHE on graphene to the present case and underline the difference. 

It is clear that, (i^p*) 2 is written in terms of the diagonal form of the Hamiltonian (|32p and 
therefore its spectrum can easily be obtained. Then, the eigenstates of H'^ take the form 

IW * 0,r,) D =( KM}'^ ) (53) 

and the zero mode is given by 

|{0},r?} D =| IL T" L )• (54) 

The corresponding eigenvalues are 

M 




ES t (m i ,r,) = ± 



D 



j^m< + -(»7-l)M(Af -1) 



i=l 



(55) 
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It can be projected on LLL that is realized by requiring that rrii = 0. This operation leads to the 
lowest energy for the excited Dirac states, such as 



Ejp(mi = 0,7?) = ±^( V - 1)M(M - 1). (56) 

Clearly, for rj = 1 we get a zero energy mode. This is an agreement with the fully occupied state v = 1 
see before fore the relativistic particles. 

To close this part, one can also report some discussions about the symmetry group that leaves 
invariance the Dirac Hamiltonian for the present case. These of course will offer a theory group 
approach to analyze the interacting relativistic particles. 

4.4 Gauge transformation 

Sometimes it is relevant to make the appropriate transformations in order to get some information 
about a system under consideration. This has an interst in sense that we can derive the integrability 
of the system from its old partner. For this, we are going to establish a concrete relation between the 
free and interacting Dirac operators. This can be done by considering an appropriate singular gauge 
transformation that captures the basic features the interacting term involved in different Hamiltonian's 
see before. It can also be seen as a kind of bosonization of the present system, for this we may refer 
to the reference [16,. The relation can be served as a good candidate to test the previous analysis 
concerning the interacting particle case. 

We start with a system of M-relativistic fermions of mass m in the presence of a external magnetic 
field. Without interaction, this system is described by the total Dirac Hamiltonian in the complex 
coordinates, such as 




E," (-2J- + fx 



This of course an immediate generalization of the Hamiltonian for one free particle. Therefore, its 
spectrum is nothing but a summation over that of one particle and the wavefunctions are tensor 
product of one single eigenfunctions. 

To explicitly determine the element that governs such transformation, we adopt the method used 
by Karabali and Sakita [16J. In doing so, let us consider two wavefunctions and establish a link 
between them. This can be done by considering a singular gauge transformation, such as 

9(z!,---,z N ) = U$(zi,---,z N ) (58) 

where ^ and $ are, respectively, totally antisymmetric and symmetric wavefunctions. This transfor- 
mation can be identified to an element of the unitarty group U(N) and should be written in terms of 
a suitable function O. This is 

U = e -i@{zx,-,z N ) _ ( 5 g) 



Clearly from (158f) we have some information about O. Indeed, it should be defined such that an 
interchange of a pair of variables gives the constraint 

exp [-i@(zx, ■ ■ ■ ,z h ■ ■ ■ ,zj, ■ ■ ■ , z N )] = - exp [-i@(z 1 , ■ ■ ■ , Zj, ■ ■ ■ , z iy ■ ■ ■ , z N )] . (60) 
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Consequentely, should be realized in such way that (|6U|) must be fulfilled. Thus, we map G in terms 
of our language as 

8(21, ■■■ ,z N ) = 7]ylmhi{zi - Zj). (61) 

i<j 

This function is not only interesting in the present case but also in the field theory. In fact, it can be 
linked to the relative angle between the position vectors the i th and j th particles. This is 

M 

Q(z u ---,z N ) = (r/-l)X>« (62) 
where 6ij is related to the Chern-Simons potential of i th particle through the relation 

a cs = (n - m E % = in - 1) E ^_Vp ( 63 ) 

., . ... ' i 1 j 

which was used to study the anyon systems. This mapping can offer another way to talk about the 
interacting particles in terms of the Dirac formalism. Also one could think to a make matrix model 
analysis of the present system by adopting the Susskind-Polychronakos approach )22i |23| . 

Now we have all ingredients to do our task. Indeed, since (|58p is a singular gauge transformation, 
the Hamiltonian for $ can be obtained from (|57p for ^f. Indeed, we show that 



tfif = U^HqU. (64) 

With this, we are now able to generate a specific potential that governs the interaction between 
particles. More precisely, is not hard to end up with the Dirac Hamiltonian for interacting particles 
in the presence of a magnetic field, which is 

„* = ±( & (" 2 i + f*)-*»£iU A 

D "^-E?'(2^ + ^,) + 2,E5„^ 

(65) 

This form shows how the gauge transformation can be employed to produce some kind of interaction 
from the free case. Moreover, the above mapping tells us that the spectrum and corresponding 
eigenfunctions of H'^ can be deduced from those of Hq. For instance, we have a similar relation 
to ([58]) between the states \mi,ni) for M-particles those are tensor products of single states (|14H17p 
and \{mi},7]) D . 

We conclude this part by noting that the involved gauge transformation can be generalized to 
higher order in terms of some interactions. This can be done by considering more that three-body 
correlations and therefore offers different integrability of some models. This will be interesting in the 
FQHE world because recall that this is based on such interactions to be understood. 



5 Composite fermion wavef unction 

To give an explanation of the filling factors beyond the Laughlin states, Jain has introduced the 
composite fermion formalism p3]. In fact, they are new kind of particles appeared in condensed 
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matter physics to provide an explanation of the behavior of particles moving in the plane when a 
strong magnetic field B is present. Particles possessing 2p, with p £ N*, flux quanta (vortices) can 
be thought of being composite fermions. One of the most important features of them is they feel 
effectively a magnetic field of the form 

B* = B± 2p$ p (66) 

where in our convention the unit flux is $o = 27r. Recalling the relation between the filling factor v 
and B, one can define a similar quantity for the field B*. This can be written as 

v* = 2vr-^. (67) 
B* ' 

It is clear that from (|66p . the factors v and v* can be linked to each other through 

v = . 68 

2pu* ±1 y ' 

This relation has been used to deal with different issues related to QHE in 2DEG. More discussions 
about the mapping (|68|) and its applications can be found in [Hj. 

As an immediate consequence of (|68p . we can map the anomalous FQHE in terms of IQHE in 
graphene. Indeed, it is easy to see 

r 2n + 1 . , 

Ap(2n + 1) ± 1 y J 



This result has been obtained in different contexts, for instance see |19j . It is obvious that for n = 0, 
we obtain Laughlin states and the anomalous IQHE can be recovered by fixing p = 0. Moreover, (|69p 
tells us that the Jain's series is quite different from the 2DEG. Therefore, it is interesting to focus on 
FQHE in graphene. 

Remembering the Laughlin potential established before, it is natural to ask about that behind 
. To reply this question, we build a wavefunction corresponding to in terms of a generalized 
potential. In the beginning, let us see how CF picture does look like for the Laughlin states, which 
are the same for both systems: semiconductors and graphene. Indeed, following Jain idea [13], one 
can write the Laughlin wavefunction as 

M / M \ 

c ug h(^ *) = n & - ex p - 4 e i^i 2 ( 7 °) 

i<j \ i / 

where xi is the wavefunction of the fully occupied state v = 1. Now it is clear that for higher Landau 
level, one may generalize (|70p to that of the form 

M ( B M \ 

^Ja?n( Z >^) =Yl ( z i ~ z j) 2P X2(2n+l)exp f -— |^| 2 J (71) 
i<j \ i J 

and X2(2n+i) is describing the composite fermions at the anomalous integer filling factor f c f . This is a 
obvious mapping between the particles generating FQHE in graphene and composite fermions behind 
unconventional integer QHE. 

At this level, one may ask about the potential of interaction that is responsible for the emergence 
of u G and therefore the model that describes such case. One way to answer this inquiry is to make a 
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straightforward generalization of the Laughlin potential (|29H30|) . Then by analogy, we can define the 
Jain potential as 

(Vj a ' ,n )i = -2i^-(lnF Jain - lnF Jain ), (W ain )i = 2i^(lnF Jain - lnF Jain ). (72) 

OZi OZi 

Clearly, for re = we recover (I29ti30p . Thus, we can write the Jain wavefunction in terms of the 
function F Ja '"(z,z) as 

^(z,z)=F^(z,z)e W . (73) 

Of course it should be described by some Hamiltonian that involves (V^ am )i and (Vg Jam )j as important 
ingredients. This requires a deep attention because of the complexity of ^'^(z, z) and therefore needs 
to be done separately. Also another important point that is one could think about a relation between 
these potentials and those of Chern-Simons in similar way to the Laughlin case. 



6 SU(N) wavefunctions 

Using different theoretical arguments, some authors suggested a possible FQHE in graphene. As we 
have seen before, this will be natural if we look at the anomalous IQHE as a product of collective 
behavior of the composite fermions instead of particles. Moreover, different wavefunctions have been 
proposed to describe the anomalous FQHE for different filling factors. Among them, we cite that 
constructed by Goerbig and Regnault [12] to solve some problems brought by other theories like for 
instance that developed in [9j. However, this wavefunction is sharing many common features with our 
early proposal [23]. On the other hand, the Goerbig and Regnault states are a direct extended version 
of those of the Laughlin [18J as well as Halperin |15| . To get more advantages, let us return back to 
the constructed SU(N) wavefunctions [13] describing the anomalous FQHE at the filling factor 

v = QiKr/qj (74) 

where is an ^ x iV matrix and qi is a vector. In fact, this included different fractions suggested 
recently for FQHE in graphene and allowed to make contact with different proposals. 

We start by reviewing the building of such wavefunctions. Indeed, to construct a general wave- 
function, we use the obtained result so far. In doing so, we consider two sectors labeled by (m) and 
(n). This is the case for instance in graphene where there are two subsystems forming a honeycomb. 
Le us define ■0( rra,n ) as a tensor product 

^(m,n) = ^(m) % ^(n) (75) 

where each ■0( m ) is given in (|17p . Assuming that the condition between matrix elements HQ,- = Kji is 
fulfilled, a natural way to construct the required wavefunction is 



N - K -K N 



nn "mri 



UT,\ TT I h-i M /M /(to) ,(m)l" mm TT I h—JM /(«)/(") /( 

I *) Jel = [[ |e 1 Mm Vfe ' %i J 1 1 [ e ni n 

m=l n=l 

N R 

n ^■■■*-+-<-vir ) • • • mn i°>- w 
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It is not hard to get the constraint 

£tot|*>Jel = 
where the total angular momenta is given by 

N 



Ltot — y ] ^ . 



(77) 



(78) 



This wavefunction is actually capturing some interesting effect that is relevant to discuss about FQHE 
in graphene. Indeed, the following term 



M 

n [«-■ 



kM m +M n ^( m ' n ),l,( m ' n ) . . . ^( m ' n ) 



k 2 



1>Mm+M„ 



Kn 



(79) 



is nothing but is an inter-layer correlation. In conclusion, our configuration could be a good ansatz 
for the ground states of FQHE in graphene. This will be clarified soon. 

At this level, we have all ingredients to do our job. Indeed, we start by defining a new complex 
variable 







(m) 



» 



,2M 



(80) 



fori = l,---,M 
for i = M + 1, • 

assuming that the particle numbers are equal, i.e. Mi = M2 = M, and recalling the antisymmetric 
Vandermonde determinant for the fully occupied state 



n 

i<3 



det ( z. 



M-j 



e h-iM z ...z M - 1 . 

U %M 



(81) 



Therefore, (j76|) can be projected on the complex plane as 

N 



Jel 



n 

N 

n[< 







M-l 



K —K N 

n 

n=l 



<■»-'» (zf 



(n) 



M-l 



■ ■ ■ f 



2M-1 

k2M 



Kn 



(82) 



It can be written in the standard form as 

K 



»- = n n 

m=l i<j 



z i z j 



nn 



An) _ (n) 
Z i Z j 



N M 

nn 

m<n ij 



(m) (n) 



A'„ 



(83) 



This exactly coincides with that constructed by Goerbig and Regnault [12] and what we have pro- 
posed [23] in terms of the matrix model and non-commutative Chern-Simons theories. The wave- 
function (|83|) is a good candidate for describing the anomalous FQHE in graphene. As it is seen 
in [13], the present wavefunction is general and has many features. In particular, it can be fixed to 
recover different filling factors like for instance those in [151 125] as well as others. Note that, the above 
wavefunctions are the planar limit of those constructed on two-sphere |13] . 

As far as the graphene is concerned, we restrict our analysis to the case N = 4. In fact, we 
give an explicit realization of the Jain potential and therefore the Hamiltonian that captures such 
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interactions. In doing so, we can use the same procedure as for the Laughlin case to write down the 
potential behind (|83p . For instance, the first term in (|83p is originated from 



l(V*%}\ Kmm = -2^(lnFi-l m - lnFi-t J, [iyi%]\ Kmm = ^(InF^ - lnF^J (84) 

and the potential corresponding to the second term can be obtained in similar way, i.e. changing K mm 
by K nn . However, the inter-layer interaction is 

{( V * e % Kmn = - 2i ^ F kL ~ toFjgj, [(Vi%\ Kmn = 2^(lnF#l - Infgj. (85) 

These can be worked as for the previous cases to get the corresponding wavefunctions. In doing so, 
we end up with the form 

i>f e r' K ™(z, z) = II F^ nm F^Fi<L *o. (86) 

m=l 

Of course this is associated to a Hamiltonian, which involves such potentials. Indeed, this can be ob- 
tained as two copies of H™ 1 amended by inter-layer correlation. It can straightforwardly be generalized 
to include other kind of interaction, in particular more that three-body type. 

The above wavefunctions has many features. Indeed, it can be linked to the composite Dirac 
fermions and derive some interesting results. To clarify this point, we return to the definition of the 
filling factor and write for any j the relation 

"> = w B < 87 » 

where the number Nb is the quantized flux. In our case, it can be written in terms of elements of the 
matrix K (I74p as 

N B = p j K j + J2Pi K H (88) 

and we have set Kjj = Kj. This offers different discussions about the above wavefunctions and its 
links to other proposals. For instance, we can establish a bridge between the filling factor (j74H and 
that coresponds to the composite Dirac fermions. For this, let us start by identifying ([69]) and (fT4"j) to 
see what are the consequences. Clearly, for n = we end up with the Laughlin series, namely 

W, = *h (89) 

which is in agreement with different investigations [111 113] . However, the interesting case is when 
p = that leads to the relation 



t^Pj 3 V2n + 1 ; 



(90) 



This offers different discussions about the ratio ^ and in particular its quantization. For this, let us 
explicitly evaluate the involved series, which can be done by recalling a similar formula. This is [26] 

> = 1 (91) 

k + 1 k\{m-k)\ m + 1 v ' 
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Therefore, by comparing both series we can identify 2n to m and more importantly the ration can be 
quantized as follows 

Pi 2n! . 

Pj ~ 2(fc + 1) k\{2n - k)\ ij ' 1 ' 

This mapping shows how the present theory can be restricted to make a possible connection with the 
composite Dirac fermions. Moreover, other theories can be recovered from our the present results like 
for instance that proposed in the reference [H]. On the other hand, some discussions can be added 
about the nature of the Hall droplet in terms of the ratio and in particular the quantization of its 
area. 



7 Conclusion 

The present work was devoted to analyze the behaviour of an interacting system in the presence of 
an external magnetic field. More precisely, we have focused on a specific case where the system is 
exactly a sheet of graphene. To do this, we have started by considering a model that captures a kind 
of interaction, which applied before at different occasions to discuss the conventional QHE. Before 
doing so, we have presented some discussions about the effective potential that corresponds to the 
Laughlin wavefunction at the filling factor v = • This was a typical example of a generic case that 
is the Jain wavefunction. 

By considering a Hamiltonian of the present system where the interaction is restricted to be two 
and three-body types, we have employed a previous technical to achieve our goal. This was done 
by making use the mapping between three Hamiltonian's: Landau, Pauli-Schrodinger and Dirac. 
This means that we have derived an appropriate Dirac operator in terms of the involved interaction. 
The mapping allowed us to get in a simple way different spectrum's entering in the game and the 
excited states as well. Using a gauge transformation, we have established a link between the free and 
interacting Dirac operators. This can be generalized to take into account of other interaction types. 

As concerning QHE in graphene, in the beginning we have extended the Laughlin potential to the 
Jain case. Indeed, we have generalized the interaction terms that are behind the series v = T^qrf to 
that for the composite fermions in graphene, i.e. v G . This was shown another way to talk about the 
Jain wavefunction for graphene and its corresponding effective potential. Of course these can be used 
to deal with other proposals related to different wavefunctions suggested for the anomalous FQHE. 

As second part, we have concentrated on the SU(N) wavefunction that captures different filling 
factors and generalizes some others. In the first step, we have given the corresponding potentials as 
an explicit example to illustrate our general procedure for Jain wavefunction. To explicitly show its 
relation to some theories, we have established its link to the composite fermion pictures. This has 
been done by identifying both of filling factors. As consequence, we have found a quantized term that 
include the density of interacting particles. This in principle gave rise to another way to interpret our 
results and made contact with others. 

The obtained results suggest to investigate and deal with other issues related to QHE in graphene. 
As example, we are now focussing on a second problem that involves another type of interaction. More 
precisely, we are wondering to understand the behaviour of the confined relativistic particles in the 
presence of an external magnetic field. 
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